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The low energy behaviour of the isotropic t-J ladder system is investigated using exact diagonaliza- 
tion techniques, specifically finding the Drude weight, the charge velocity and the compressibility. 
By applying the ideas of Luttinger liquid theory, we determine the correlation exponent Kp which 
defines the behaviour of the long range correlations in the system. The boundary to phase separation 
is determined and a phase diagram is presented. 

At low electron density, a Tomonaga-Luttinger-like phase is stabilized whilst at higher electron 
densities a gapped phase with power law pairing correlations is stabilized: A large region of this 
gapped phase is found to exhibit dominant superconducting correlations. 

PACS numbers: 74.72.-h, 71.27.+a, 71.55.-i 
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I. INTRODUCTION 

Over the last few years, the behaviour of strongly cor- 
related electrons confined to coupled chains has received 
widespread attention; the reasons for this are numer- 
ous. Firstly, the behaviour of electrons in one-dimension, 
under t-J or Hubbard type interactions, is now rela- 
tively well understood and described generally by the 
term Luttinger Liquid (LL). The coupling of two such 
Luttinger liquids as in the ladder geometry provides an 
interesting first step towards the challenge of describ- 
ing the behaviour in two-dimensional systems. A sec- 
ond reason for interest in these systems lies in the un- 
usual nature of the ground state in the undoped sys- 
tem, namely spin liquid behaviour with a finite gap in 
the spin excitation spectrumlll; this behaviour is in con- 
trast to the gapless behaviour of a single chain. The 
evolution of the spin-gapped state on doping has obvious 
relevance to gapped superconducting behaviour. In acU 
dition, compounds such as (1^02)^2070 and SrCu20^ 
are believed to be well described by a lattice of coupled 
chains. Very recently, experiments on Lai-xSrxCu02.^ 
have provided insight into the doping of coupled chain 
systems. Whilst there is considerable literature on many 
aspects of the t-J ladder behaviour, a complete picture 
is still far from being realised: Our aim in this paper is 
to clarify some of the behaviour of the t-J ladder sys- 
tem by drawing on some of the ideas used in describing 
the strictly one-dimensional systems (i.e. LL theory). In 
conjunction with results from several other techniques, 
we will then present a speculative phase diagram. 

The t-J Hamiltonian on the 2x L ladder is defined as, 
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where most notations are standard. /3 (=1,2) labels the 
two legs of the ladder (oriented along the a;-axis) while j 
is a rung index (j=l,...,L). We shall concentrate on the 
isotropic case where the intra-ladder (along x) couplings 
J and t are equal to the inter-ladder (along y) couplings 
J' and t'. 

At half filling the hamiltonian reduces to the Heisen- 
berg model|-and the behaviour is generally relatively well 
understoodtl. A simple interpretation is given by con- 
sidering the strong coupling limit ( J = 0) in which the 
ground state consists of a singlet on each rung with a 
spin gap J') which corresponds to forming a triplet 
on one of the rungs. With the introduction of intra-chain 
coupling J, the triplets can propogate and form a coher- 
ent band thereby reducing the spin gap. In the isotropic 
case, the spin gap remains (~ 0.5J). 

The evolution of this spin gapped state on doping is 
perhaps one of the most interesting aspects of the ladder, 
behavioum. Recent work on this hole-doped phasd3iJ'E3 
has indicated a finite spin gap, a single gapless charge 
mode, hole pairing and possible dominant superconduct- 
ing correlations. In this paper we shall discuss some new 
independent results which provide a more complete de- 
scription of not only this gapped phase but the whole 
region of parameter space. A possible phase diagram for 
the isotropic t-J ladder as^a function of J/t and doping 
has been proposed recentlycl and we use these ideas in our 
analysis: Away from half filling the spin-gapped phase 
is stabilized up to J/t ~ 2.1 where the system phase 
separateai3. As the system is doped further, a phase 
with a single gapless spin and a single gapless charge 
mode is found and, as we shall explain, this behaviour 
is like that of the one-dimensional Tomonaga-Luttinger- 
Liquid system. As in the gapped phase, phase separation 
occurs for large J/t. At very small electron densities, an 
electron paired phase exists. Note that although in the 
t-J ladder there are four possible zero-momentum gap- 
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less modes (two spin and two charge), throughout the 
phase diagram only one gapless charge mode and either 
zero or one gapless spin modes are observed; this then 
allows an almost identical treatment to the strictly one- 
dimensional case. 

In section |l| we discuss briefly the Luttingcr Liquid 
theory used to describe strictly one-dimensional systems 
and also the possible application of this theory to coupled 
chains; In section [J] we present our numerical calcula- 



the variation of the ground state energy Eq with particle 
density n, i.e. the inverse compressibility 



tions, and finally in section [V we apply the Luttinger 
Liquid theory to our results and present a phase diagram 
for the system. 



II. LUTTINGER LIQUID BEHAVIOUR: 
TOMONAGA-LUTTINGER AND 
LUTHER-EMERY PHASES 

In dealing with strictly one-dimensional interacting 
fermion systems, cme can in general make use of con- 
formal field theoryEj and bosonizatiorJlj which allow a 
determination of the decay exponents of the various cor- 
relation functions to be determined from the low energy 
behaviour of the model. The general idea is that one- 
dimensional interacting fermion systems can be mapped 
onto the Fermi-gas model and the corresponding 'g-ology' 
weak coupling theoryll3. This Fermi gas model scales to 
two different regimes, namely the Tomonaga-Luttinger 
(TL) fixed point and the Luther-Emery (LE) line which 
are relevant for repulsive {gi > 0) and attractive {gi < 0) 
backscattering matrix elements respectively; as we shall 
explain, the important difference between these two uni- 
versality classes lies in the spin degrees of freedom. The 
low lying excitations of the Fermi gas model are collec- 
tive spin or charge density oscillations, which propagate 
with different velocities, giving rise to spin-charge sepa- 
ration and power law behaviour of the correlation func- 
tions. In the TL phase, both a gapless spin and a gapless 
charge mode are exhibited; In contrast, whilst exhibiting 
a gapless charge mode, the LE phase has a gap to spin 
excitations. 

Conformal field theory relates the properties of a fi- 
nite system (such as the compressibility and the Drude 
weight) with the correlation exponents; One coefficient 
(Kp ) determines the exponents of all the power law de- 
cays, (and similarly the singularity of the momentum 
distribution function close to fc/). Hence if a particu- 
lar model scales to the TL (or LE) universality class, one 
can infer the dominant correlations from the low energy 
behaviour of the system which can be deduced from much 
smaller system sizes than would be required to calculate 
the correlation lengths directly. 

The relationships between the correlation exponent 
(Kp) and the low energy behaviour of the model are given 
below for a system of size N, and length L ( we have cho- 
sen to give general equations such that N=L for a chain, 
N=2L for the ladder geometry) . Firstly, the ratio of the 
charge velocity Up to the coefficient Kp is proportional to 
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The coefficient Kp is also related to the Drude weight 
(To; this Drude weight is the weight of the zero frequency 
(dc) peak in the conductivity a^} and may be obtained 
by considering the curvature of the ground state energy 
level as a function of threaded flux, $ 



CTo = 'i.UpKp = 
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We can also determine the charge velocity by considering 
the dispersion of the energy spectrum 



[E^p - Eo)/{2t/L) 



(4) 



where Eip is the lowest lying charge mode to the ground 
state {Eq) with neighbouring k value. 

These equations provide us with 3 independent con- 
ditions on Kp and Up which can be used to check the 
consistency of the Luttinger liquid relations. Also a cal- 
culation of the parameter Kp is relatively straightforward 
and this parameter then determines the exponent coeffi- 
cients of all the correlation functions. 

As explained previously, the essential difference be- 
tween the TL and LE fixed points lies in the spin de- 
grees of freedom: the LE region gapped whilst the TL 
gapless. The correlation exponents of the two different 
cases are summarised in table I (taken from referenced). 
We have omitted the logajrithmic corrections and these 
are detailed in referencdlJ. We emphasise that for LL 
systems, the correlation functions show either power law 
or exponential decay, with interaction-dependent pow- 
ers determined by one coefficient Kp. Also we see that 
for Kp < 1 (spin or charge) density waves at 2k f are 
enhanced and diverge, whereas for Kp > 1 pairing fluc- 
tuations dominate. 

Whilst in strictly one-dimensional systems there is a 
single gapless charge mode, the theory could equally well 
be applied in a case where more than one gapless mode 
existed if the excitations were decoupled in the low energy 
regime. In such a case, each degree of freedom would have 
an associated operator algebra, i.e. an associated Kp. 

Specifically considering the problem of coupled chains, 
we note that at present little is known and there is much 
interest into how to connect the quasi-onc-dimensional 
results to the strictljp-pne-dimensional case. A recent 
calculation by SchulzEj has considered the coupling of 
two Luttinger liquids by a small interchain hopping us- 
ing a bosonization technique. Interestingly, in the pres- 
ence of both forward and backward scattering terms, 
the calculation predicts a gap in all the magnetic ex- 
citations and a gapless charge mode (as observed for 
the hole-doped region of the t-J ladder). This gapped 
phase however exhibits somewhat different correlations 
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to the strictly one-dimensional LE phase described above. 
Firstly the CDW^r and SDW.,r correlations decay expo- 
nentially along the chains (we use the notation of Schulz 
where or tt indicate the oscillations are in or out of 
phase between the two chains respectively, i.e. 'bonding' 
or 'antibonding'). A divergent density-density response, 

decaying as ~ cos (^2(fcy -I- kj)rjr~^'' , exists in analogy 

with the 4k f oscillations of a single chain {k'j and kj refer 
to the Fermi points of the bonding and antibonding quasi 
particle branches respectively); we note that the coeffi- 
cient Kp differs by a factor of two to that used by Schulz 
since we have chosen to define ctq and k as the Drude 
weight and compressibility per site rather than per rung. 
The superconducting correlations (cross chain pairing) 
decay as r"^/^*" and exhibit a 'd'-like character. Hence 
for this gapped phase, we would expect dominant super- 
conducting correlations for Kp > 1. In agreement with 
these findings, different approaches bAt-|Troyer et a£3, 
NagaosaO and and Balents and Fisherlla predict similar 
behaviour in the spin-gapped region: i.e. 'Luther-Emery- 
like' in the sense that there exist two order parameters 
(analogous to the on-site pairing and 2k f CDW in the LE 
class) whose exponents obey a reciprocal relation: these 
correspond to the pair field correlations and to a four 
fermion operator {nB{r)nB{0)) , where is the density 
of 'bosonic' hole pairs bound on a rung (analogous to the 
Akf oscillations of a single chain {i.e.2k^+2kj)). In table 
II we summarize the correlation exponents predicted for 
this spin gapped phase in the ladder geometry. 

In the following section we give details of calculations 
using these ideas to characterize the behaviour of the t-J 
ladder, finding the parameter Kp and hence the dominant 
correlation functions. We consider various electron den- 
sities and various ratios of J/t to build up a speculative 
phase diagram. 

III. NUMERICAL CALCULATIONS 

Since we require information concerning the low energy 
properties of the model, the dominant technique we have 
employed is that of exact diagonalization of finite sys- 
tems, specifically 2x5 and 2 x 10 double chain rings. 
Exact diagonalization techniques are particularly well 
adapted to the investigation of low energy modes since 
implementation of various quantum numbers is straight- 
forward; the various excitation modes can be obtained 
by calculating the ground state energy in each symmetry 
sector. The low energy modes of the system are charac- 
terized firstly by their spin: singlet and triplet excitations 
correspond to charge and spin modes respectively. It is 
also useful to consider the parity of the states under a 
reflection in the symmetry axis of the ladder along the 
direction of the chains: Even {R^ = 1) or odd {Rx = —1) 
excitations corresponding to bonding (B) or anti-bopding 
(A) modes respectively (0 and tt as used by Schulzlij) . Fi- 
nally the dispersion relation of each mode is determined 



by the momentum kx=2'Kn/ L. 

In order to ensure that the antiferromagnetic correla- 
tions are not frustrated when one goes around each chain, 
we have chosen the electron number to be always a multi- 
ple of four; our results then concern electron densities 0.4 
and 0.8 for both system sizes and in addition 0.2 and 0.6 
for the larger system size. The absolute ground state is 
given by the boundary conditions that form a closed shell 
in the non-interacting Fermi sea (obtained by turning off 
the interaction, J) and these are used in the calculations 
of Up and k, specifically anti-periodic boundary condi- 
tions for n < 0.5 and periodic boundary conditions for 
n > 0.5. 



A. Drude Weight and Anomalous Flux Quantization 

The first calculation we present concerns the Drude 
weight, defined by equation 0. The numerical technique 
involves threading the double chain ring with a flux $ 
and studying the functional form of the ground state en- 
ergy with respect to the threaded flux, namely i?o($). 
In general Eo{^) consists of a series of parabola, cor- 
responding to the curves of the individual many body 
states En{^): This envelope exhibits a periodicity of one, 
where we have chosen to measure the flux in units of the 
flux quantum $0 = hc/e. Note that the function i?o(*&) 
also gives a quantitative value of the sujprfluid density 
Dg, which is in general different from ctq^- The Drude 
weight corresponds to curvature of a single ground-state 
many-body energy level whilst the superfluid density cor- 
responds to the curvature of the envelope of the individ- 
ual many-body states as a function of flux. However, 
since the flux {(jjc) at which another many body energy 
level crosses the zero-flux ground state energy level varies 
as 0c ~ {he/ e)L^^'^ (where d is the dimension )li3, in one- 
dimension 4>c is independent of L; there are only a finite 
number of energy level crossings in the the thermody- 
namic limit and (Tq and Ds are equal (up to a factor of 
2tt). 

In addition to the Drude weight and the superfluid 
density, the function Eo{^) also yields information re- 
garding the phenomenon of anomalous flux quantization; 
this has been explained in a previous publicationLl so we 
mention it only briefly here. Whilst in general the ground 
state envelopje-|i?o(<I>) exhibits a periodicity of one, By- 
ers and Yang23 have shown that in the thermodynamic 
limit Eq(<^) exhibits local minima at quantized values of 
flux, the separation of which is 1 /n where n is the sum 
of the charges in the basic group. Hence, for a paired 
superconducting state we would expect minima in Eq(<^) 
at intervals of 1/2. These minima are related to the ex- 
istence of supercurrents which are trapped in metastable 
states correspondin|g-|to the flux minima and are thus un- 
able to decay awaynil. It should be mentioned that this 
anomalous flux quantization (AFQ) is an indication of 
pairing and is not in itself sufficient to imply a supercon- 
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ducting state. 

Numerically the application of a flux through the dou- 
ble chain ring is achieved by modifying the kinetic term 
of the hamiltonian such that 



(5) 



where $ is the flux through the ring measured in units 
of $0- Hence the application of a flux is numerically 
equivalent to a change in the boundary conditions of the 
problem; $ = representing periodic and $ = 1/2 rep- 
resenting anti-periodic boundary conditions. In the ther- 
modynamic limit, (To must be iudependent of the phase 
introduced at the ring boundaryE3 and therefore we con- 
sider the whole of the envelope Eo{^) as a function of 
flux (in general consisting of several parabola). 

Choosing the parameters n — 0.8 and J/t ~ 0.5, we 
show in figure |l|a(b) all the possible spin and charge 
modes of the 2 x 5 (2 x 10) system, for all possible mo- 
menta, as a function of applied flux. In the case of the 
larger system, we show the full spectrum for $ < 0.25 in 
order to simplifM the diagram (this work has been previ- 
ously publishedu). For both system sizes the minimum 
energy is formed by charge (spin zero) bonding modes; 
the excited modes with different quantum numbers move 
further from the ground state as the system size is in- 
creased (a result we have checked by finite size scaling) 
and hence will not interfere with Eo{^). The existence of 
minima at intervals of half a flux quantum (i.e. anoma- 
lous flux quantization) clearly indicates the existence of 
pairing). 

The envelope L[Eo{^)~ Eq{^ — 0)] has been extracted 
and is shown in figure 0a(b) along with equivalent plots 
from other regions of the phase diagram: Figure j^a shows 
the data for J/t =1.0 and electron densities of 0.4 and 
0.8, whilst figure |^ shows the data for an electron den- 
sity of 0.8 for ratios J/t = 0.5 and 4.0. Apart from the 
curve with n = 0.8, J/t = 4.0 (which appears to scale to 
a flat function, consistent with a phase separated state), 
all the data appears to show only small finite size effects. 
Anomalous flux quantization is observed for the larger 
electron density (for the lower values of J/t) indicating 
pairing, and hence consistent with a superfluid state in 
this region. 

In order to determine the Drude weight, we simply 
calculate the average value of the curvature of L[Eo{^) — 
Eo{^ — 0)] over all <&; a quadratic curve was fitted to 
each portion. In figure ^a we plot the Drude weight as 
a function of J/t for electron densities 0.2,0.4,0.6 and 
0.8 for the 2 x 10 system and electron densities 0.4 and 
0.8 for the 2x5 system. The curves are plotted up to 
a maximum in J/t which is determined by the value at 
which the system phase separates (see compressibility). 
In figure ^ we plot the Drude weight as a function of 
electron density for various values of the ratio J/t. 

There are several features of the resulting behaviour 
we should mention: Note firstly that finite size effects are 
relatively small with the 2x5 results close to those of the 



2 X 10 results. The Drude weight increases as the electron 
density is increased from zero, until it reaches one-quarter 
filling, then decreases with increasing electron density. 
Also as we would expect for a spin charge separated state, 
the Drude weight is effectively independent of J. 



B. Charge Velocity 

The second quantity we have calculated is the charge 
velocity, defined by equation ^. Considering the charge 
bonding modes (the lowest lying charge modes), the en- 
ergy difference between the ground state energy level and 
the energy level with neighbouring momentum, Akx = 
^ was calculated. As an example, we show in figure ^ 
the charge bonding modes for the case n = 0.4 J/t = 2.0, 
indicating with solid lines the specific energy levels whose 
energy difference gives the charge velocity. We note that 
the gap (and therefore Up) is approximately constant as 
a function of fiux. For consistency with future data how- 
ever, we have calculated the charge velocity at the par- 
ticular flux which gives the absolute ground state, i.e. 
$ = 0.5 for n < 0.5 and <i> = for n > 0.5; The results 
are shown in figure ^ as a function of J/t for various 
system sizes and various electron densities (the results 
for n = 0.6 are not included since the numerics present 
some difficulty close to one-quarter filling). Note again 
that finite size effects are relatively small. 



C. Compressibility 

The next quantity we calculate is that of the compress- 
ibility, defined by equation |^. The finite size equaivalent 
is given by 
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Eo{n + An) + Eo{n - An) - 2Eo{n) 
(An)2 



(6) 



where Eo^n) is the ground state energy of the finite sys- 
tem of ladder length L with an electron density n. As 
for the calculation of the charge velocity, the boundary 
conditions have been chosen to give the absolute ground 
state. An represents the finite change in electron density, 
0.2 and 0.4 for the 2 x 10 and 2x5 systems respectively. 
In Fig. ^ we show the results of the calculation of 1/n^K 
for electron densities corresponding to 0.2,0.4,0.6 and 0.8 
for the 2 X 10 ladder and we also plot the result of the 
2x5 system for an electron density of 0.4. 

In addition to the determination of the specific val- 
ues of the inverse compressibility, the boundary to phase 
separation may also be determined from these results. 
It is well known that at sufficiently large values of J/t, 
a system will undergo a separation into two phases; a 
hole-rich and an electron rich phase. This effect arises 
principally to minimize the number of broken antiferro- 
magnetic bonds in the system. Phase separation occurs 
when the compressibility diverges (i.e. a 'liquid' to a 
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'solid' phase) and hence the inverse compressibility van- 
ishes. Whilst there remain some small finite size effects, 
the general form of the phase separation line is readily ob- 
served from figure ^. As electron density is increased, the 
value of J/t at which phase separation occurs decreases 
(although electron densities of 0.6 and 0.8 are both indi- 
cate phase separation close to J/t ~ 2.1). -These results 
agree well with those of Tsunetsugu et alE3 who used a 
similar technique but varied both system size and elec- 
tron density simultaneously. The phase separation curve 
can be extrapolated to all electron densities and will be 
shown in the predicted phase diagram, figure |[ 

IV. LUTTINGER LIQUID PARAMETERS FOR 
THE LADDER 

In order to explore the validity of the Luttinger liquid 
relations in our problem, we consider the ratio do/Trn^KWp 
which equals unity for a Luttinger liquid. The results of 
the numerical calculation of this quantity are shown in 
figure H for various electron densities. 

At low electron densities (i.e. for the cases n = 0.2 and 
n = 0.4), previous work has suggested that a TL phase is 
stabilized and the results of this ratio show good agree- 
ment with the predicted value of unity; for the case of 
n = 0.4 an increase in system size shows the ratio scaling 
towards unity. For the hole-doped spin-gapped region 
(n = 0.8), where the behaviour is less well understood, 
the system phase separates at a much lower value of J/t 
and hence the curve drops to zero at J/t ~ 2.1. Before 
this phase separation, the data is not inconsistent with 
the 'new' LE-like behaviour which may be described by 
equations (finite size effects are largest for high elec- 
tron density and low J/t). Since the ratio is close to 
unity, it appears to confirm the earlier justification for 
using this one-dimensional theory, i.e. only a single gap- 
less charge mode is observed. 

With the values of the compressibility and the Drude 
weight, we can obtain an estimate of the coefficient Kp 
using Kp = l/2y^7Tn'^KaQ. This behaviour of Kp for the 
different electron densities of the 2 x 10 ladder is shown 
in Fig. and we have also plotted the results for a 2 x 5 
ladder with electron density 0.4. Note that as J/t is in- 
creased, Kp increases (for all electron densities) becom- 
ing infinite at phase separation as the compressibility 
diverges. A-sLmilar calculation has been performed by 
Troyer et aOilj for a specific electron density of 0.857, 
i.e. two holes on a 2 x 7 ladder and the results are con- 
sistent in this region. 

Before analysing the behaviour further, specifically the 
importance of Kp, we present a speculative phase dia- 
gram of the isotropic t-J ladder as a function of J/ 1 and 
electron density. This phase diagram is shown in figure 
P and we discuss briefly the various regions. 

At larger values of J/t, the system phase separates, 
and to estimate the value of J/t at which this occurs. 



we show the data points at which the inverse compress- 
ibility vanishes in the 2 x 10 system (see figure ||). On 
doping away from half filling, the spin gap region per- 
sists and a phase exhibiting one gapless charge mode is 
stabilized. On further doping a gapless phase is stabi- 
lized: the schematic boundary between these two phases 
is shown as a dot-dashed line. As for both the one- and 
two- dimensional t-J casesEj, a gas of electron pairs is 
formed at low electron densities above a critical value of 
the ratio J/t; other 2p-particle (p>l) boundstates could 
also become stable at larger J/t in this region. Again 
the boundary to this paired phase is dot-dashed and is 
schematic. 

From the data in figure ^ we have plotted contours of 
constant Kp in order to allow a determination of the dom- 
inant correlation functions. Whilst the parameter Kp is 
continuous for a-^articular value of J/t as the electron 
density is variedE3, at some region between n — 0.4 and 
n — 0.8 a gap opens in the spin excitation spectrum and 
the correlation functions change their form discontinu- 
ously, scaling to a different fixed point as explained in sec- 
tion ^ A 'jump' in the exponent of the superconducting 
correlations occurs with the SCd charge exponent chang- 
ing from 1 -I- 1/Kp to l/Kp-, the 2k f CDW correlations 
jump from power law behaviour (exponent 1 -I- Kp) to 
exponential decay. In addition, in the gapped high den- 
sity state we would expect conjugate 'four fermion 4fc^' 
CDW correlations with exponent Kp. 

For both the high density 'gapped' phase and the low 
density TL phase we expect different correlation func- 
tions to dominate either side of the contour Kp = 1; 
in both cases superconducting correlations dominate for 
Kp > 1. However, in the gapped hole-doped phase, this 
region is much larger and in contraaLto other models such 
as the one-dimensional t-J modelo, does not just exist 
as a precursor to phase separation. A physical picture of 
the behaviour is of the holes pairing up on the rungs, the 
spins in singlets, and the dominant correlation functions 
are then associated with the movement of the hole pairs. 
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FIG. 1. Energy as a function of flux (in units of $o = hc/e) 
for a) the 2x5 and b) the 2 x 10 ladders with J/t = 0.5 and 
n = 0.8. We show all possible momenta for various quantum 
numbers: For the charge modes, the solid lines correspond to 
bonding and the dotted lines to anti-bonding, while for the 
spin modes the dashed lines correspond to bonding and the 
dot-dashed lines to anti-bonding. For the larger system size, 
we give only the charge bonding and the lowest lying spin 
antibonding mode in full to simplify the diagram. 



FIG. 2. L[Eoi$) - So($ = 0)] where L is the length of the 
ladder and Eo{$) is the ground state energy with an applied 
flux $. The dashed lines correspond to 2 x 5, the solid lines 
to 2 X 10. a) shows the results for n = 0.4 and n = 0.8 both 
with J/t = 1.0, while b) shows J/t = 0.5 and J/t = 4.0 both 
with n — 0.8. 



FIG. 4. (a) Energy of the charge bonding modes as a func- 
tion of flux for the 2 x 10 system with n = OA, J/t = 2.0. The 
solid lines indicate the ground state and the excited states 
used to calculate the charge velocity, (b) The charge velocity 
Up as a function of ratio J/t for various electron densities and 
system sizes. The dotted lines are a guide to the eye. 



FIG. 5. l/n^ function of ratio J/t for various electron 

densities and system sizes. The dotted lines are a guide to the 
eye. 



FIG. 6. The ratio ivn^Ku^p/ao as a function of J/t for vari- 
ous electron densities and system sizes. The dotted lines are 
a guide to the eye. 



FIG. 7. Kp as a function of ratio J/t for various electron 
densities and system sizes. The dotted lines are a guide to 
the eye. 



FIG. 8. Speculative phase diagram of the t-J ladder as a 
function of J/t and electron density n. The circles and crosses 
represent results from the 2 x 10 system and the dotted lines 
are guides to the eye. The dot-dashed lines separating the 
different phases are estimated. 



TABLE I. Correlation exponents for the Tomonaga Lut- 
tinger and Luther Emery cases. SS and TS indicate singlet 
and triplet superconductivity (pairing) correlations respec- 
tively. 



Correlation 


TL 


LE 


2k f SDW 


l + Kp 


exponential 


2k f CDW 


l + Kp 


Kp 


SS 


1 + l/Kp 


l/Kp 


TS 


1 + l/Kp 


exponential 


4fc/ CDW 


4Kp 


4A-p 



TABLE II. Correlation exponents for the spin gapped 
phase in the ladder geometry. The n indicates that the cor- 
relations along the two chains are out of phase. 



Correlation 


gapped spin 


2k f SDWvr 


exponential 


2k f CDW^ 


exponential 


SCd 


yKp 


4fc/ CDW 


Kp 



FIG. 3. The Drude weight ao a) as a function of ratio J/t 
for various electron densities and system sizes, and b) as a 
function of electron density for various J/t. The smooth lines 
are a guide to the eye. 
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